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Abstract 

We report on a fully self-consistent determination of a phase transition to 
a superconducting state in a conserving approximation. The transition tem- 
perature calculated for a two-dimensional Hubbard model with an attractive 
interaction in the fluctuation exchange approximation agrees with quantum 
Monte Carlo calculations. The temperature dependences of the superfluid den- 
sity and of the specific heat near T c indicate that the phase transition in this 
model of coupled collective degrees of freedom and electronic degrees of free- 
dom is consistent with neither mean-field theory, Gaussian fluctuations about 
a mean field order parameter, nor unbinding Kosterlitz-Thouless vortices. 
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The quasi two-dimensional nature of copper-oxide superconductors has focused 
considerable attention on the role of fluctuations associated directly with supercon- 
ductivity 0. In a strictly two-dimensional superconductor 0, thermal fluctuations 
of the order parameter preclude a finite temperature phase transition to a state with 
long-range order |||, f|. A Kosterlitz-Thouless transition may occur between two 'dis- 
ordered' phases, with the low temperature phase characterized by a finite superfluid 
density and correlation functions that decay algebraically with distance ||. Should 
fluctuations be important for the copper-oxide superconductors, their effects would 
be most apparent in the normal state where weak three-dimensional coupling does not 
substantially influence the essentially two-dimensional character of electronic proper- 
ties. Fluctuations of the superconducting order parameter may explain the pseudogap 
in the normal-state single-particle excitation spectrum as well as the anomalous tem- 
perature dependences observed for transport and thermodynamic properties in the 
normal state pi 0]. 
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Figure 1: Internal energy (diamonds) and free energy (squares) as a function of 
the threaded flux calculated in the FEA for a 2D Hubbard model on a 32 x 32 
lattice with U = —4 and n = 0.75 for a temperature above (open) and below (filled) 
the superconducting transition. The energy and free energy are measured with 
respect to their zero-flux values. Below the superfluid transition, an energy barrier 
at 4>b/4>o = 0.25 appears for both E and F for the single period shown. 

To explore the role of Cooper-pair fluctuations, we focus on the evolution of the 
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superconducting state with decreasing temperature in the fluctuation exchange ap- 
proximation (FEA) 0, a conserving approximation || beyond mean field theory. 
The FEA and the Hubbard Hamiltonian || with an attractive interaction consti- 
tute a minimal model of coupled quasiparticles and s-wave Cooper-pair fluctuations. 
We report the first observation of a phase transition in the full FEA. The FEA su- 
perconducting phase transition occurs at a temperature in excellent agreement with 
quantum Monte Carlo (QMC) calculations. 

We probe superconductivity by calculating the internal energy E and free energy 
F as a function of magnetic flux <pB threading the periodic lattice. Byers and Yang ]10], 



TT| argued that off-diagonal long range order makes E(4>b) for Cooper pairs periodic 
in 4> B with period 4>o/2 = hc/2e. Scalapino, White, and Zhang |12| subsequently 
showed that the superfluid density determines the curvature of F(<J)b) near <p B — 0, 

F(<p B )^F(0)+ 1 -D s (T)(^y + ..., (1) 

where D s (T)/2ti 2 is the superfluid density in units where the hopping is t = 1 (see 
below). Similarly, the energy E((J)b) satisfies 

E(<PB)*F(0) + ±iD s -T^](^f + ..., (2) 

which has been used to argue for a KT transition from Monte Carlo calculations 
on small lattices. Explicit quantum Monte Carlo calculations of E(<Pb) as a func- 
tion of temperature for the two-dimensional attractive Hubbard model show a low- 



temperature superfluid state [y| [L4| consistent with the onset of algebraic pair cor- 
relations [jn| |nj . 

Fig. [p illustrates one of our central results: E{(pB) and F((p B ) are flat at high 
temperature, but at low temperatures become periodic with period <j) /2 and barrier 
at o /4, indicating the formation of a low-temperature superconducting state fTTH - 



We consider a Hubbard model on an L x L lattice with periodic boundary condi- 
tions. In an applied magnetic flux <pB, the Hubbard Hamiltonian is, 

H{(j> B ) = ~t [exp[ 2 ^ B ] 4 c r+ - CT + h.c] 

r,a Vo 1 ^ 

-t [ c Lcr+y, CT + h.c.} + nr,in T ,i, (3) 

r,cr r,cr 

where —t is the nearest-neighbor hopping matrix element and U (< 0) is an on-site 
(attractive) interaction. 

We calculate the grand thermodynamic potential Q, from the (self-consistent) 
self-energy S and propagator G using the Luttinger-Ward formula fL8f , 



Q (T, n, <j) B ) = -2 Tr [ EG + ln(— G" 1 + E) ] + $[G], (4) 
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where 3>[G] generates the skeleton diagram expansion for the self-energy, 
which, together with Dyson's equation, provides closed equations for G and E. 



Complete expressions for &fea[G] and T,fea[G\ appear elsewhere \T§\. The FEA 



for £ contains exchanged spin, density, and Cooper-pair fluctuations [^(J. The con 



tribution from Cooper-pair fluctuations is £ pp (r, r) = T pp (r, r) G{— r, — r), where 

r pp (q,a; m ) = pp 6 

is an effective non-local and retarded interaction generated by repeated interactions of 
s-wave pairs, and Xpp( r > r ) — G(r, r) G(r, r) is the particle-particle bubble. We do not 
explicitly break gauge symmetry, and hence no anomalous one-particle propagators 
enter our calculations, but off-diagonal order (of arbitrary spin or orbital symmetry) 
can still appear spontaneously in our two-particle propagators, because the functional 
$ is explicitly of infinite order in U. 

The thermodynamic description of a superfluid transition requires reliable calcula- 
tions of small differences in internal energies and free energies for different These 
are calculated from G, S, and Q, 

E(<Pb) = ^Tr [{2e k + E}G], (7) 
F(4> B ) = Q + /m, (8) 

where e k is the bare dispersion relation, and as in Eq. [|, all quantities on the 
right-hand side of Eqs. (0) and (|), including e^, are functions of 4>b- We use a 
Fourier-transform-based parallel algorithm to solve self-consistently the FEA equa- 



tions I^T], ^] . Contributions to convolution sums from high-frequency tails, essential 
for accurate calculations of thermodynamic properties, are included analytically. 

Calculating D s (T) from Eq. ([D using the fully self-consistent propagators and free 
energy is equivalent to obtaining D S (T) from the full conserving FEA current-current 
correlation function; hence our observation of superconductivity is not inconsistent 
with previous work that did not see evidence for superconductivity in an approximate 
FEA response function [23]. 



The energy shown in Fig. [T] was calculated for a 32 x 32 lattice with U = —4 
and a density n = 0.75 P7[ . Fig. |2] shows the temperature evolution of the energy 
barrier A^E = E{(p /A) — E(0) for small lattices in comparison with QMC results 
|13|| and mean field theory (MFT). Our results for n = 1.0 show no evidence for 
superconductivity, consistent with QMC calculations and theoretical expectations. 
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The negative values for A^E in both the QMC and FEA and the downturn at low-T 
in the FEA are apparently finite-size effects. For n = 0.75, the results from FEA and 
QMC agree surprisingly well, but the increase with lattice size in the maximum of 
A^E, interpreted as a signature of a KT transition in the QMC results, is significantly 
smaller in the FEA. 

n=0.75 n=1.0 
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Figure 2: Temperature-dependent energy barriers A^E, for U — — 4 calcu- 
lated in the FEA (filled diamonds) shown in comparison with the quantum 
Monte Carlo data of Ref. 11 (open circles) for lattices of the same size, and 
A^E calculated in mean-field theory (open squares). The FEA shows a finite 
barrier (indicating superfluidity) for n = 0.75 (left column) and a transition 
temperature well below the mean-field theory value of ~ 0.7 and in good 
agreement with the transition temperature ~ 0.1 obtained in quantum Monte 
Carlo calculations for lattices of the same size. Like the quantum Monte Carlo 
calculations, we find no evidence for superfluidity in the FEA for n = 1. 

The QMC calculations focussed on A^E, the barrier at 0o/4, because these cal- 
culations have significant noise and the response is largest for this 4>b ■ In our calcu- 
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lations, numerical noise is much smaller and we are able to use smaller values of the 
flux to estimate the superfluid density and its temperature derivative from Eqs. (U) 
and (|2|). The top panel of Fig. (|^) shows our estimates of D S {T) and dD s {T) / dT 
for a 32 x 32 lattice obtained with a flux of 4>b/4>o — 0.1; errors due to higher-order 
terms in Eqs. ([I]) and (0) are expected to be less than 2.5%. 
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Figure 3: (top) The estimated superfluid density D S (T) (squares), and 
dD s (T)/dT (circles) as obtained from the flux dependence of the free and 
internal energies for a 32 x 32 Hubbard lattice with U = —4 and n = 0.75. 
The slowly increasing D S (T) with decreasing temperature does not suggest the 
form D s oc Q(Txt — T) expected for the vortex-driven transition of Kosterlitz 
and Thouless. (bottom) The onset of superfluidity is evidently associated with 
a near-instability of the particle-particle fluctuation T-matrix as signaled by 
the approach of — Ux PP (ci = 0, u m = 0) to unity. 



The behavior of D S (T) shown in Fig. (^) differs from expectations for both a 
mean-field transition (D s oc T c — T below T c in the thermodynamic limit) and for a 
KT transition (D s is discontinuous at T c in the thermodynamic limit). On a 64 x 64 
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lattice, the rise in D S (T) occurs at slightly lower temperature, but D S (T) is otherwise 
very similar to the result shown for a 32 x 32 lattice. 

The bottom panel in Fig. [3] shows that the superfluid transition occurs as T pp 
approaches an instability, i.e. as U\pp approaches —1. The T-matrix is not the self- 
consistent pair-susceptibility for the FEA and a near instability does not in general 
signal a phase transition. In this case, extrapolation of — Ux PP (<i = 0,c<j m = 0) for 
T ~ 0.25 to unity yields a temperature in good agreement with that for the onset of 
energy-barrier formation. 
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Figure 4: Specific heat as a function of temperature for a 32 x 32 site 2D 
Hubbard lattice (diamonds) with U = —4 and n = 0.75. The temperature of 
the peak in the specific heat correlates with the transition temperature. Also 
shown is the specific heat for a 64 x 64 lattice (open circles) for the same 
density and interaction strength. 



A prominent peak appears in the specific heat |24| as shown in Fig. f|. The 
peak is not consistent with MFT, which predicts a jump in the specific heat, nor 
with Gaussian fluctuations about a mean-field order parameter. In the latter case 
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fluctuations round the corners of the jump in reasonably narrow critical region. The 
specific heat is also not consistent with expectations for a KT transition, which shows 
a peak above T c from the entropy associated with unbinding topological vortices. 
Comparison of the specific heat and entropy shown for 32 x 32 and 64 x 64 lattices 
indicates that finite-size effects in the results shown for these quantities are small. 

The superconducting transition we observe in the FEA is apparently neither a 
simple mean-field transition nor a KT transition, unless the FEA results for D S (T) 
contain unusual finite-size effects that die very slowly with increasing lattice size. 
Our tentative explanation is that the anomalous temperature dependence we observe 
in these quantities results from the more complicated interplay between collective 
degrees of freedom (order parameter fluctuations) and electronic degrees of freedom 
in the FEA than in either the BCS mean-field theory (which describes the change 
in single-electron energies but ignores order-parameter fluctuations) or descriptions 
based on Ginzburg-Landau functionals (which treat fluctuations of the order param- 
eter, but ignore their effects on electronic quasiparticles). Although we do not mean 
to suggest that the details of the results for D S (T) or C(T) from the FEA are correct, 
we do think that the FEA provides a valuable example of how calculated properties of 
a superconductor can be modified by mutual feedback between order-parameter fluc- 
tuations and quasiparticle properties when the feedback is included self-consistently. 
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